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Preview 



This is a review devoted to some results of Algebraic Programming 
(Computer Algebra - CA) used in treating several problems of general 
relativity, based mainly on already published articles (see the refer- 
ences) . 

The first chapter presents procedures in REDUCE language using 
EXCALC package for algebraic programming of the Hamiltonian for- 
mulation of general relativity (ADM formalism). The procedures cal- 
culate the dynamic and the constraint equations and, in addition, we 
have extended the obtained procedures in order to perform a complete 
ADM rcductional procedure. Several versions of the procedures were 
realized for the canonical treatment of pure gravity, gravity in interac- 
tion with material fields, infiationary models (based on a scalar field 
non-minimally coupled with gravity) and theories with higher order 
lagrangians. 

The second chapter is devoted to the same problem as the first one, 
but now we use Maple V platform. In this purpose we present Maple 
procedures using the GrTensorll package adapted for the study of the 
canonical version of the general relativity based on the ADM formalism. 

Last chapter presents CA procedures and routines applied to the 
Dirac field on curved spacetimes. The main part of the procedures is 
devoted to the construction of Pauli and Dirac matrices algebra on an 
anholonomic orthonormal reference frame. Then we shall present how 
these algebraic programming sequences are used in the study of the 
Dirac equation on curved spacetimes and noninertial reference frames. 
A comparative review of such procedures obtained for the two men- 
tioned CA platforms (REDUCE + EXCALC and MAPLE + GRTen- 
sorll) is presented. Applications for the calculus of Dirac equation 
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CONTENTS 



on specific examples of spacetimes witli or witliout torsion and for tlie 
study of tlie non-relativistic approximation of tlie Dirac field are pointed 
out, including the search for inertial effects. 



Chapter 1 



REDUCE programming and 
the Hamiltonian version of 
general relativity 

1.1 Introduction 

The Hamiltonian formalism of the general relativity is based on the 
(3+l)-dimensional split of space-time - ([|^, 0, [0. For instance, the 
Hamiltonian formulation and quantization of some homogeneous and 
inhomogeneous cosmological models needs the complete expressions 
of the super-Hamiltonian Ti and super-momentum IH} (see below) and 
the dynamic equations in term of the canonical variables. For this 
purpose, concerning the great volume of calculations and the division of 
the calculus in distinct steps it is possible to create algebraic procedures 
which transpose the specific manipulations of the canonical formulation 
of general relativity in computer language. 

There are some early results in this direction [Q, obtained in old 
versions of the REDUCE language, without EXCALC package. Tra- 
ditionally, the programs calculate only the dynamic equations and the 
constraint equations after introducing the canonical variables and their 
canonically conjugate momenta and do not perform the complete re- 
ductional formalism in order to point out the true dynamic content of 
the treated model. 
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In the present chapter we shall summarize recent results obtained 
processing some space-time models with several new procedures realized 
with EXCALC package (in REDUCE). In fact, here we have a review of 
a series of articles - ||^) containing the details of our procedures 
and the obtained results for some concrete space-time models. 

We have extended the programs in order to realize a complete re- 
ductional procedure (solving the constraint equations, changing of vari- 
ables, reduction of dynamic variables, etc.) for some space-time models. 

We shall present also the version of our procedures for the Hamil- 
tonian treatment of some inflationary models (based on a scalar field 
non- minimally coupled with gravity). 

1.2 The local form of the canonical for- 
malism of gravity 

Here we shall use the specific notations for the ADM formalism 
for example latin indices will run from 1 to 3 and greek indices from 
to 3. The starting point of the canonical formulation of the general 
relativity is the (3+l)-dimensional split of the space-time produced by 
the split of the metric tensor : 



^ ^ ^ 9oo ^ ^ 9oj 

(4), 



( NkN'' - N^ N, 



\ ^ ^ 9io ^ ^ 9ij 



J 



XI] 



\ Ni gij 



where Qij is the riemannian metric tensor of the three-dimensional 
spacelike hypersurfaces at t = const, which realize the spacetime folia- 
tion. Here N is the " lapse" function and A^* are the components of the 
"shift" vector §. 

The Einstein vacuum field equations now are (denoting by "•" the 
time derivatives) : 

Qij = 2Ng-'/'[7c,j - ]^g^,^l\] + N,,^ + N^/, (1.2) 



= -Ng'/^[R^ - ^g^m] + ^Ng-^l^g^^i-K^'TXki - \{^\)^] 
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+ [n'^N"']/^ - N'/^t:^^ - ,mT^^' (1.3) 

where vr*-' are the components of the momenta canonically conjugate to 
the gij's. 

In the above formulas we denoted by "/" the three-dimensional 
covariant derivative defined with gij using the components of the three- 
dimensional connection F*^ and Rij is the three-dimensional Ricci ten- 
sor. 

The initial data on the t = const hypersurface are not independent 
because they must satisfy the constraint equations, which complete the 
Einstein equations 

n = -^{R + g-^[^-{n\y - vr*%,,]} = (1.4) 

H' = -27r'^/j = (1.5) 

where Ti. is the super-Hamiltonian, the super-momentum and g is 
the determinant of the three-dimensional metric tensor gij. The action 
functional in Hamiltonian form for a vacuum space -time can thus be 
written as (]l[,[g) 

S = j dt j {n'^gij -NH- N.H'p^uj^uj^ (1.6) 

where the tu^'s are the basis one-forms. 

In Dirac's version of the Hamiltonian formalism p| the constraints 
equations ( p..4| ) and (|1.5|) are not solved and no coordinate condition is 
imposed; The Dirac Hamiltonian is then given by 

H = J {NH + NiH'pWuj^ (1.7) 



Thus the dynamic equations ( |1.2|) and (|1.3| ) are obtained by differ- 
entiating H with respect to the canonical conjugate pair of variables 



(TT'^S-fcr, 



When interaction with a scalar field x is considered, the free grav- 
itational action must be completed with several new terms. The la- 
grangian of the massless scalar field is [|l^ Lsc = — (— '■^''fi')^^^ g^''X,^iX,v 
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and defining the scalar field momentum as tt^ = {dLgd dx,t) we must 
add to the total gravitational action : 



Ssc = j dtjin^x-Nij^ + y^g'^x,Xj)-Nii7T^g''x,j))uj'uj^'^S (1.8) 

It is obvious from the above formula that the terms in brackets 
must be added to the gravitational super-Hamiltonian and to the super- 
momentum. Thus we can obtain the new form of the constraint equa- 
tions and Dirac Hamiltonian. With the last one it can obtain the 
dynamic equations for the scalar field variables. 

When interaction with a perfect fluid (consisting of pressureless 
dust) is present, to the total gravitational action we can add : 

Sd= fdtf (vr^ - Nn^ \l + g'^^i(P,j] ' - iV,(?^^0,,-)^'^'^' (1-9) 



where is the scalar field describing the dust. 

In the last few years, inflationary models became very popular not 
only in cosmology but also in the theory of general relativity and grav- 
itation. The majority of the mentioned models are based on the inter- 
action of gravity with a scalar field [|I5| - [[l^ , , ||38[ . We shall treat 
a model based on the action functional : 

S = j d^x^^g (^$2 - 1 + A$^^ (1.10) 

where the scalar field $ is non-minimally coupled with the gravity. 
As an inflationary model the above action functional determines the 
behavior of the gravitational constant as being G = l/lQn^"^ through 
the dynamics of the field $. 

In order to obtain the canonical version of the field theory described 
by the action functional ( |1.10| ) we must redefine the components of the 



three-dimensional momentum canonical conjugate to the 

n'^ = U:>K - K'A (1.11) 



where Kij = i^i/j + ^j/i ~ 9ij) is the extrinsic curvature of the 
three-dimensional hypersurfaces. With K we have denoted the trace of 
Kij, namely K = K^i. 
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Thus, after an appropriate calculation presented in complete details 
in we obtain the complete canonical version of our inflationary model 

S = Jdt (fx (rr'^gij + 7r$<i> - NH - NiH') (1.12) 

where the new super-Hamiltonian Ti and super-momentum Ti^ (and the 
new constraint equations) are: 

^^g'^<l>,<l>, + 47^(7'^$$,,, - ^A<^' = (1.13) 

n' = -2n'^/j - Ag'^^-^^jTT + = (1.14) 

The dynamical equations are obtained in the usual way. For exam- 
ple, after defining the Dirac Hamiltonian as : 

Hd = J dt d^x {NH + NiTt) (1.15) 



we can set : 



. _5Hd _ 6Hd 



^TT*-?' 6gij 
As a result we have : 

2 1 

9ij = Niij + Nj/i + —N^'^TXij - —N^'^Txgij- 

.^N^-\^g,, - ^-^^.pNug'^'^gi^ (1.16) 
V 9 



V " 
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.^iV<l>'*$'^' - -^g'^NA<!>^ + (a^^<I>'J' + A^^$'*) (1.17) 



17 

+ yv--- - - 2 

Of course, in a similar manner we can obtain the dynamic equations 
for the scalar field $ and his conjugate momenta 7r$. We do not give 
here the mentioned equations because, in the algebraic procedures we 
shall present in the next paragraph we have obtained these equations 
directly from the Dirac Hamiltonian Hd by a variational derivative (a 



facility of EXCALC package [0]). The complete form of the gij^s and 
TT^-^'s are necessary as a consequence of the fact that in practice we use 
the components of the above tensors as functions of the true canonical 
variables introduced in every specific case of interest. 



1.3 The REDUCE computer procedures 

We shall briefly present the steps of calculus necessary for developing 
the canonical program outlined above. The steps will be the same if 
the calculus is realized manually or with the computer, using specific 
EXCALC procedures. Several steps are the same if we have pure grav- 
ity, gravity coupled with matter fields or inflationary models. Thus we 
shall present the next steps together with examples of computer lines 
from the pure gravity topic. Here are the steps : 

1) introduction of the basis vectors and coframe forms (with a usual 
metric statement in EXCALC [^) together with the calculus of the 
structure coefficients of the Cartan basis considered; it is realized, in 
the same time the introduction of the metric gtj components, because, 
in EXCALC package the assignment of the coframe forms defines, in 
the same time the metric components; 

2) introduction of the momentum vr*-^ defining the canonical conju- 
gates of the metric components and their dependence on the spatial 
variables defining a symmetric 2-form, named PIU; here is introduced 
a program line computing the term vr*-' ^^ from eq. (|1.6|) and verifying 
if the dynamic part of the action functional is in canonical form : 



act := piu(j ,k)*@(g(-j ,-k) ,t) ; 
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if the answer is negative this step must be run again, after a new defi- 
nition of the 71^^ components; 

3) evaluates the Ricci tensor Rij and the Ricci scalar R of the metric 
above defined by direct computation from a RIEMANNCONX form 
named OM (a riemannian connection form), declared in the previous 
lines of the program 

ricci (-j ,-m) := ee(-in)_ I (ee(-k)_ I (d om(k,-j) + 

oni(k,-p)"om(p,-j))) ; 

scricci := ricci(j,-j); 

where ee{—k) are the basis vectors defined together with the coframe; 

4) calculates the super-Hamiltonian Ti. : 

haO : = -sqrt (detm ! *) * (scric+d/detm ! *) * 

((l/2)*( (piu(j,-j))*(piu(k,-k)) ) 
-piu(j ,k)*piu(-k,-j) ) ); 



from equation (|1.4|) or from equation ( |1.13|) in our inflationary model 

5) calculates the super-momentum W : 

pform ha(j)=0; 
ha(j):= 

-2*((ee(-k)_|d piu(j ,k) ) + (ee(-k)_ I oni( j , -p) ) *piu(p,k) 
+ (ee(-k)_|oni(k,-p))*piu(j ,p) 
-(ee(-k)_|oni(p,-p))*piu(j ,k)) ; 

with equation ( |1.5| ) or eq. ( |1.14| ) for the inflationary model; 

6) at this stage it is already possible to calculate the temporal deriva- 
tives for the gij's components (eq. ( |1.2| ) or ( |1.16D ) : 

pform n=0,derge(j ,k)=0; tvector ni(k); 
fdomain n=n(r,t) ,ni=ni(r,t) ; 
derge(-j ,-k) := 2*n* (1/sqrt (detm! *) ) * 



12 



CHAPTER 1. REDUCE PROGRAMMING AND 



( piu(-j,-k)-(l/2)*g(-j,-k)*(piu(p,-p)) ) + 
(ee(-k)_|d ni(-j)) - (ee (-k)_ I om(p, -j ) ) *ni (-p) + 
(ee(-j)_|d ni(-k)) - (ee (-j )_ I om(l , -k) ) *ni (-1) ; 

(where n and ni are the "lapse-shift" functions) and then the evolution 
equations of the canonical variables defined in step 1; 

7) calculates the second set of temporal derivatives for 7r^\ with eq. 
( p..3| ) (or ( |1.17| )) and the dynamic equations for momenta defined in 



step 2; 

All the above mentioned steps are naturally transposed in EXCALC 
lines, as we can see from the above examples. Special program lines 
are devoted to the calculus of the new terms appearing in the theory, 
when interaction of gravity with a matter field is considered. It is 
calculated the new Dirac Hamiltonian (obtained adding new terms to 
the super-Hamiltonian and the super-momentum specific for the matter 
field). Using Hamilton- type equations, the dynamic equations of the 
scalar field variables (time derivatives of x, and vr^, tt^) are obtained 
together with the new terms added to the older dynamic equations. 

We have realized also file-sequences in order to perform the reduc- 
tional ADM formalism for the canonical theory presented above. There 
is not an unique method to realize this purpose. The procedure must 
be adapted to the specific space-time model . The point is to guess a 
canonical transformation (or a sequence of canonical transformations) 
which, after solving the constraint equations and imposing specific co- 
ordinate conditions, generates an action with one or two dynamic vari- 
ables and a Hamiltonian generating the time behavior of the system 
without constraints. These attempts to perform a complete reduc- 
tional formalism in each processed model can be a good way to realize 
the connection with numerical relativity. We investigate the possibility 
to generate FORTRAN lines for numerical solving of the differential 
equations obtained after the canonical program outlined here is per- 
formed. 

When inflationary models are considered the above program lines 
are adapted accordingly to the new dynamic and constraint equations 
(as we already mentioned - see also P]) and, in addition we have new 
steps devoted to the dynamical treatment of the scalar field. We shall 
use a different method as a consequence of the fact that the scalar field 
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has a unique component : 

8) it is calculated the Dirac Hamiltonian from the super-Hamiltonian 
and supermomentum obtained in steps 4 and 5 with formula ( 1.15| ). 



9) the dynamical equations for $ and 7r$ are obtained using a facility 
of EXCALC package, namely the variational derivative : 

derphi := vardf (had,piphi) ; 
derpiphi := - vardf (had, phi) ; 

where "derphi" and "derpiphi" denote the temporal derivatives of the 
scalar field $ and of his momentum 7r$ respectively. With "had" we 
have denoted the Dirac Hamiltonian already mentioned. 



1.4 About the concrete results 

In the above cited articles (|^ - [0) we have presented several results 
obtained with space-time models processed by the procedures presented 
here. 

In order to verify the accuracy of our programs we approached some 
well known models from the literature 0,0, [1^ - [1^. 

We have reobtained all the results reported in 0]. Some of these 
space-time models were processed in a generalized form in interaction 
with a massless scalar field. Another well known case, the canonical 
treatment of cylindrical gravitational waves |jl8| (also coupled with a 
scalar field) was one of the preferred testing models for us. The coinci- 
dence of our results with the well-known results from [|18| was a good 



sign for approaching other models, more sophisticated. Among these 
models is also the space-time with three-dimensional subspace which 
we shall present here. We give the form of the spatial metric tensor and 
of the momenta canonically conjugate to the Qi/s, the constraints (the 
super-hamiltonian and the super-momenta) as well as the expressions 
of the dynamic equations for all the canonical variables. 

At the end of the list of specific results we have presented the skele- 
ton of the reductional procedure we have applied here : changing of 
variables, canonical transformations, rescaling of variables, coordinate 
conditions etc., and finally the reduced hamiltonian. 
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Space-time model with a three sub-space in form of a 3- 
torus here depicted as a 2-tori fibring of 



9ij 



( e2("'-/) 
e-2(-+/) 
V e^" 



where a — a{z,t), w — w{z,t), f — f{t) and < x,y, z < 1, and the 
coframe is — dx, w"^ — dy, — dz. 





V 











_ie2{-+/)(7r^ + 7r^) o 







2^ "a 



The derivatives with respect to z variable being denoted with "/" 
we have 



a = N^a' + N^' + ^c^-^-^TYtt^ + Je^^-^iVTr/ ; x = iV^x' + ^e^^^'NTr^ 



nf = N'n'f + N^'TTf - -Ne'f-'^{nl + 27r„7r/ + tt^) 
-4e-2/-»(a'iV' - iV" - Nw'^) - ^e-2^-"7V(e^^7rJ - 4x") 
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o 

+ \e-^-'^N{e'f7rl + 4x'') 



- In order to perform a reductional ADM formalism we have some 
different possibilities, but the most cases lead to nonlocal theories. A 
natural and very simple version is suggested by Misner, (in - but not 
used in the cited article) is to choose / as a time variable (the transverse 
cross-sectional area of the universe is proportional to e^-^). Taking as 
coordinate variable w, solving the constraint equations in terms of vr^ 
and 7Tf and imposing the ADM coordinate conditions T = f = t and 
w = z we can obtain the action functional as : 



5 = y dtdz{nad + n^x ~ Hred) 



where the reduced time dependent hamiltonian Hred is 



27r„ 



Now, in the final lines, we shall mention only the results obtained 
after processing the computer procedures outlined in the previous sec- 
tions for some space-time models, considered as initial data for the 
inflationary universe described by ( 1.10 ). The results reported in |0 



are only in a brute form without any regarding of the physical sig- 
nificance. Spherically-symmetric model (usual in the dynamic study of 
inflationary universes), Bianchi I/Kantowski-Sachs models and also the 
above model with T^ subspace were presented in details in H . 
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Chapter 2 



Using Maple in the study of 
the canonical formalism of 
general relativity 



2.1 Introduction 

The use of computer facilities cam be an important tool for teaching 
general relativity. We have experienced several packages of procedures, 
(in REDUCE + EXCALC for algebraic programming and in Mathe- 
matica for graphic visualizations) which fulfill this purpose. In this 
chapter we shall present some new procedures in MapleV using GrTen- 
sorll package ([^) adapted for the canonical version of the general 



relativity (in the so called ADM formalism based on the 3+1 split of 
spacetime). This formalism is widely used in the last years as a ma- 
jor tool in numerical relativity for calculating violent processes as, for 
example the head-on collisions of black holes, massive stars or other 
astrophysical objects. Thus we used these computer procedures in the 
process of teaching the canonical formalism as an introductory part 
of a series of lectures on numerical relativity for graduated students. 
We shall use the same notations and conventions already presented in 
the first chapter. Obviously we used the programs in REDUCE pre- 
sented there for producing our new procedures for Maple + GrTensorll 
package, but because there are many specific features we shall present 
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here in some detail these procedures in the next section of the chapter. 
The last section of the chapter is dedicated to the conclusions pointed 
out by running the Maple procedures presented here and some future 
prospectives on their usage toward the numerical relativity. 



2.2 Maple + GrTensorll procedures 

Here we shall describe briefly the structure and the main features of the 
Maple procedures for the canonical formalism of the general relativity as 
described in the previous chapter. Two major parts of the programs can 
be detected : one before introducing the metric of the spacetime used 
(consisting in several definitions of tensor objects which are common 
to all spacetimes) and the second one, having line-commands specific 
to each version. 

The first part of the program starts after initalisation of the GrTen- 
sorll package (grtw();) and has mainly the next lines : 

> grdefCtr := pi{"i i}'); 

> grdef ( 'haO : =-sqrt (detg) * (Ricciscalar+ 

(l/detg)*((l/2)*(tr)-2-pi{i j>*pi{ "i }))'); 

> grdef ('ha{ }:=-2*(pi{ "i "j ;j>-pi{ ~i "j >*Chr{ p j "p » ' ) ; 

> grdef ('derge{ i j }: =2*N(x,t) * (detg) " (-1/2) *(pi{ i j } - 

(l/2)*g{ i j}*tr)+Ni{ i ;j } + Ni{ j ;i }'); 

> grdef ('Ndd{ "m j > : = Nd{ "m ; j > ' ) ; 

> grdef ('bum{ "i ~j ~ni>:=pi{ "i ~j }*Ni{ ~in }'); 

> grdef ('bla{ ~i ~j }:=bum{ "i "j "m ;m }'); 

> grdef ( 'derpi{ ~i "j }:= 

-N(x,t)*(detg) " (1/2)*(R{ "i "j }-(l/2)*g{ ~i "j }*Ricciscalar)+ 
(l/2)*N(x,t)*(detg)-(-l/2)*g{ -i -j >*(pi{ "k "1 }*pi{ k 1 }- 
(l/2)*(tr)~2)-2*N(x,t)*(detg)-(-l/2)*(pi{ "i "m >*pi{ ~j m >- 
(l/2)*pi{ -i -j }*tr)+ (detg)- (1/2) *(Ndd{ "i "j }-g{ "i "j }* 

Ndd{ "m m >) + bla{ ~i "j > - Ni{ ~i ;m }*pi{ "m ~j }- 

Ni-[ ~j ;m >*pi{ "m "i }') ; 

Here haO and ha{ ' i } represents the superhamiltonian and the 
supermomentum as defined in eqs. (|1.4| ) and (|1.5| ) respectively and tr is 



the trace of momentum tensor density vr*-' - which will be defined in the 
next lines of the program. Here N(x,t)represents the lapse function A^. 
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Also, derge{ i j } represents the time derivatives of the components 
of the metric tensor, as defined in eq. ( |1.2| ) and derpi{ " i " j } the time 
derivatives of the components of the momentum tensor tt*-' as defined 
in eq. (|1.3|). 

The next hne of the program is a specific GrTensorll command for 
loading the spacetime metric. Here Maple loads a file (previously gener- 
ated) for introducing the components of the metric tensor as functions 
of the coordinates. We also reproduced here the output of the Maple 
session showing the metric structure of the spacetime we introduced. 



> qloadCCyl.din') ; 



Default spacetime = Cyl_din 
For the Cyl_din spacetime: 
Coordinates 
x(up) 

a 

X = [x, y, z] 
Line element 
2 2 
ds = exp(gamma(x, t) - psi(x, t)) d x 

2 2 2 

+ R(x, t) exp(-psi(x, t)) d y + exp(psi(x, t)) d z 

As is obvious we introduced above the metric for a spacetime with 
cylindrical symmetry, an example we used for teaching purposes being 
a well known example in the literature ([|l^). In natural output this 
metric has the form : 



9ij 



j eT-^ 

R^e-^ I (2.1) 
V 



in cylindrical coordinates x,y,z with x G [0, oo), y G [0,27r), z G 
(— oo, +oo) where R, ip and 7 are functions of x and t only. 

After the metric of the spacetime is established the next sequence of 
the programm just introduce the components of the momentum tensor 
vr*-^ as 
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> grdef('Nd{ " m > := [dif f (N(x,t) ,x) , 0, 0]'); 

> grdef('Ni{ "i } := [Nl(x,t), N2(x,t), N3(x,t)]'); 

> grdef ('vil{"i}: = [pig(x,t)*exp(psi(x,t)-gaiiima(x,t)) ,0,0] ') ; 

> grdef ('vi3{~i> 

:=[0,0,exp(-psi(x,t))*(pig(x,t)+(l/2)*R(x,t)*pir(x,t)+ 

pip(x,t))] ') ; 

> grdef ('vi2{"i}: = [0, (2*R(x,t))"(-l)*pir(x,t)*exp(psi(x,t)) ,0] ') ; 

> grdef ('pi{ "i "j } := 

vil{ "i >*kdelta{"j $x}+vi2{ "i }*kdelta{ "j$y }+ 
vi3{ >*kdelta{~j $z}'); 

> grcalc(pi(up,up)) ; 

> grdisplay (pi (up , up) ) ; 

Here Ni{ " i } represents the shift vector A^* and the other objects 
(Nd, vil, vi2 and vi3) represent intermediate vectors defined in order 
to introduce the momenum pi{ " i " j } having the form : 



/ TT^e'^-^ \ 

^nne^ 



(2.2) 



In the program we denoted Tr-y, ttr and tt^ with pig, pir and pip, 
respectively. The momentum components are introduced in order that 

the dynamic part of the action of the theory be in canonical form, that 
is : gijn"^^ = 7r^7 + n^ip + tvrR. The next lines of the programm check 
if this condition is fullfiled : 

> grdef ('del{ i }: = [diff (grcomponent(g(dn,dn) , [x,x] ) ,t) ,0,0] ') 

> grdef ('de2{ i >: = [0,diff (grconiponent(g(dii,dn) , [y,y] ) ,t) ,0] ') 

> grdef ('de3{ i } := [0,0,dif f (grcomponent (g(dn,dn) , [z,z] ) ,t)] ' ) 

> grdef ('ddgt({ i j >:= 

del{ i }*kdelta{j $x>+de2{ i }*kdelta{ j$y }+ 
de3{ i }*kdelta{ j $z>'); 

> grcalc (ddgt (dn , dn) ) ; 

> grdef ('act: =pi{ "i "j >*ddgt{ i j }'); 

> grcalc(act); gralter(act,siiiiplify) ; grdisplay (act) ; 

By inspecting this last output from the Maple worksheet, the user 
can decide if it is necessary to redifine the components of the momentum 
tensor or to go further. Here the components of the momentum tensor 
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were calculated by hand but, of course a more experienced user can try 
to introduce here a sequence of commands for automatic calculation of 
the momentum tensor components using the above condition, through 
an intensive use of solve Maple command. 

Now comes the must important part of the routine, dedicated to 
calculations of different objects previously defined : 

> grcalc(haO) ; gralter(haO, simplify) ; 

> grdisplay (haO) ; 

> grcalc(ha(up)) ; gralter(ha(up) , simplify) ; 

> grdisplay (ha(up) ) ; 

> grcalc(derge(dii,dii)) ; gralter(derge(dn,dii) , simplify) ; 

> grdisplay(derge(dn,dn)) ; 

> dl : =exp(-psi (x, t) ) *gr component (derge (dn,dn) , [z ,z] )+exp(psi (x,t)- 

gamma(x,t) ) *grcomponent (derge (dn,dn) , [x,x] ) ; 

> simplify(dl) ; 

> d2:=(l/ (2*R(x,t)))*exp(psi(x,t))*grcomponent(derge(dn,dii) , [y,y] )+ 

(l/2)*R(x,t)*exp(-psi(x,t))*grcomponent(derge(dn,dii) , [z,z] ) ; 

> simplify (d2) ; 

> d3 :=exp(-psi(x,t))*gr component (derge (dn.dn) , [z,z] ) ; 

> simplify (d3) ; 

> grcalc(derpi(up,up)) ; gralter (derpi (up, up) .simplify) ; 

> grdisplay(derpi(up,up)) ; 

> fl := exp(gamma(x,t)-psi (x,t) ) *grcomponent (derpi (up.up) , [x,x] )- 

pig(x,t)*(d3-dl); 

> simplif y(f 1) ; 

> f2:= 2*R(x,t)*exp(-psi(x,t))*grcomponent(derpi(up,up) , [y,y] )+ 

(l/R(x,t))*d2*pir(x,t)-pir(x,t)*d3; 

> simplify(f2) ; 

> f3 := exp(psi(x,t))*grcomponent(derpi(up,up) , [z,z] )+d3*(pig(x,t)+ 

(l/2)*R(x,t)*pir(x,t)+pip(x,t))-f l-(l/2)*R(x,t)*f2- 
(l/2)*pir(x,t)*d2; 

> simplify(f3) ; 

This is a simple series of alternation of grcalc, gr alter and grdis- 
play commands for obtainig the superhamiltonian, supermomentum 
and the dynamic equations for the theory, dl ... d3 and fl ... f3 

arc the time derivatives of the dynamic variables, 7, R, ip and tt^, ttji, 
TT^ respectively. Denoting with "/" the derivatives with respect to r we 
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display here the resuhs for the example used above (cylindrical gravi- 
tational waves) : 

= e^{2R" - i?Y + \{i^'?R - n.TTR + ^(vr^)') = 

j = N'^-f' + 2N'^' - e^NiTR ; R = R' - Nn^ 

K 

TT^ = iVV; + N^'tt^ - e^{R!N' + \r!^'N - \^''^RN + ^iYvr^vri? - ^N^l) 

2 A 2 4x1 

jr^ = N^n'ji + N^'ttb, + - 2N" - 2N'ip' + -j'x(j'N - x(j"N - -0'^ + —Ntt 

2 2R 

jr^ = N^tt'^ + N^'tt^ + e^iRN'il:' - R"N + ^NR'j' + R'll^'N - ^-f'lp'NR 

These are the well-known results reported in (|]18|) or ([^]). 

One of the important goals of the canonical formalism of the general 
relativity (which constitutes the "kernel" of the ADM formalism) is the 
reductional formalism. Here we obtain the true dynamical status of the 
theory, by reducing the number of the variables through solving the 
constraint equations. This formalism is applicable only to a restricted 
number of space-time models, one of them being the above cylindrical 
gravitational waves model. Unfortunately only a specific strategy can 
be used in every model. Thus the next lines of our program must 
be rewritten specifically in every case. Here, for teaching purposes we 
present our example of cylindrical gravitational wave space-time model. 
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Of course we enccourage the student to apply his own strategy for other 
examples he dares to calculate. 

In our example of cylindrical gravitational waves, the reductional 
strategy as described in ([0) starts with the usual rescahng of H and 
W to n and W by 

n = e^n ; N = e^N ; Ti} = e'-^lH} ; N^ = e^-^N^ 
wich produce the next sequence of Maple+GrTensorll commands : 

> grdef ( ' ahaO : =sqrt (exp (gamma (x,t)-psi (x,t) ) ) *haO ' ) ; 

> grdef ('aha-[ "j } := exp(gainma(x,t)-psi(x,t))*ha{ ~j }'); 

> grdef ( ' an: =sqrt (expCpsi (x,t)-gainma(x, t) ) ) *n(x, t) ' ) ; 

> grdef ('ani{ "i } := exp(psi(x,t)-gamma(x,t))*ni{ ~i }'); 

The canonical transformation to the new variables, including Kuchar's 
"extrinsic time", defined by : 

T = T(oo) + l\-'K,)dr , Ht = -y + [In {{Hf - (T')')]' 

J CO 



R = R , Iln = 7rn + [\n{^j±^)]' 

are introduced with : 

> pig(x,t) :=-diff (T(x,t) ,x) ; 

> pir(x,t) :=piR(x,t) - dif f (ln( (dif f (R(x,t) ,x)+dif f (T(x,t) ,x) )/ 

(diff (R(x,t) ,x)-diff (T(x,t) ,x))) ,x) ; 

and specific substitutions in the dynamic objects of the theory : 

> grmapChaO, subs , dif f (ganima(x,t) ,x) = 

diff( ln( (diff (R(x,t) ,x))-2- (dif f (T(x,t) ,x) ) "2 ) ,x)-piT(x,t) , 'x') ; 

> grcalc(haO) ; gralter(haO,simplify) ; 

> grdi splay (haO) ; 

> grmap(ha(up) , subs , dif f (gamma(x,t) ,x)=dif f ( ln( (dif f (R(x,t) ,x) ) "2- 

(diff (T(x,t) ,x))-2 ) ,x)-piT(x,t) , 'x') ; 

> gralter(ha(up) , simplify) ; 
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> grdisplay (ha(up) ) ; 

> grcalc(ahaO) ; 

> grmapCahaO, subs , diff (gamma (x,t) ,x)=diff ( ln( (dif f (R(x,t) ,x) ) "2- 

(diff (T(x,t) ,x))-2 ) ,x)-piT(x,t) , 'x') ; 

> gralter(ahaO, simplify, sqrt) ; 

> grdisplay (ahaO) ; 

> grcalc(aha(up)) ; 

> grmap(aha(up) , subs , dif f (ganma(x,t) ,x)=dif f ( ln( (dif f (R(x,t) ,x) ) "2- 

(diff (T(x,t) ,x))~2 ) ,x)-piT(x,t) , 'x') ; 

> gralter(aha(up) , simplify) ; 

> grdisplay (aha(up) ) ; 

> grmap(act, subs , dif f (gamma (x,t) ,x)=diff ( ln( (dif f (R(x,t) ,x) ) "2- 

(diff (T(x,t) ,x))-2 ) ,x)-piT(x,t) , 'x') ; 

> grcalc(act); grdisplay (act) ; 

Thus the action yields (modulo divergences) : 

/oo roo , _ _ _ _ 

dt / driUrT + HrR + n^^j + tt^x - NH - N^Ti}) 
-oo ^0 

where : 

= t'Ht + rHIr + v'tt^ + ^T^x 

Solving the constraint equations ?i = and 

li} = for and Hr and imposing the coordinate conditions 
T = t and R = r wc obtain finally : 

/-|-oo r+oo ]_ 

-oo ^0 ^ 

from the next sequence of programm lines : 

> R(x,t):=x; T(x,t):=t; grdisplay (ahaO) ; 

> solve (grcomponent (ahaO) ,piT (x , t) ) ; 

> piT(x,t):= -l/2*(x"2*diff (psi(x,t) ,x)"2+pip(x,t)"2)/x; 

> eval(piR(x,t)) ; 
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> piR(x,t) :=-diff (psi(x,t) ,x)*pip(x,t) ; piR(x,t) ; 

> grdisplay (ahaO) ; grdisplay(aha(up) ) ; 

> piT(x,t) ; 

2 /d \2 2 

x I — psi(x, t) I + pip(x, t) 
\dx / 

- 1/2 

x 

> piR(x,t) ; 

/d \ 
-| — psi(x, t) I pip(x, t) 
\dx / 

> grcalc(act); grdisplay (act) ; 

For the Cyl_din spacetime: 

act 

/d \ 
act = I — psi(x, t) I pip(x, t) 
\dt / 

> grdef ('Action:=act+piT(x,t)*diff (T(x,t) ,t)+piR(x,t)*diff (R(x,t) ,t) ') ; 

> grcalc (Action) ;gralter(Action, factor, normal, sort, expand) ; 

> grdisplay(Action) ; 

For the Cyl_din spacetime: 
Action 



/d \2 /d \ 

Action = - 1/2 X I — psi(x, t) I + I — psi(x, t) I pip(x, t) 
\dx / \dt / 
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2 

pip(x, t) 
- 1/2 

X 

> grdef ('Ham:=piT(x,t)*diff (T(x,t) ,t)+piR(x,t)*dif f (R(x,t) ,t) ') ; 

> grcalc(Ham); gralter (Ham, expand) ; 

> grdisplay (Ham) ; 

For the Cyl_din spacetime: 
Ham 

2 

/d \2 pip(x, t) 

Ham = - 1/2 X |~ psi(x, t) I - 1/2 

\dx / X 



2.3 Conclusions. Further improuvements 

We used the programms presented above in the computer room with the 
students from the graduate course on Numerical Relativity. The main 
purpose was to introduce faster the elements of the canonical version 
of relativity with the declared objective to skip the long and not very 
straitforward hand calculations necessary to process an entire example 
of spacetime model. We encouraged the students to try to modify the 
procedures in order to compute new examples. 

The major conclusion is that this method is indeed usefuU for an 
attractive and fast teaching of the methods involved in the ADM for- 
malism. On the other hand we can use and modify these programs for 
obtaining the equations necessary for the numerical relativity. In fact 
we intend to expand our Maple worksheets for the case of axisymmet- 
ric model (used in the numerical treatement of the head-on collision of 
black-holes). Of course, for numerical solving of the dynamic equations 
obtained here we need more improuvements of the codes for paralel 
computing and more sophisticated numerical methods. But this will 
be the object of another article. 



Chapter 3 

Gravity, torsion, Dirac field 
and computer algebra 



3.1 Introduction 



In a series of recent published articles |^-||33[ we have presented some 
routines and their applications written in REDUCE+EXCALC alge- 
braic manipulation language for doing calculations in Dirac theory on 
curved spacetimes. Including the Dirac fields in gravitation theory re- 
quires lengthy (or cumbersome) calculations which appropriately could 
be solved by computer algebra methods. Initially our main purpose 
was to develop a complete algebraic programming package for this pur- 
pose using only the REDUCE -|- EXCALC platform. Partially this 
program was completed and the main results and applications were re- 
ported in our above cited articles |3^-[^. But we are aware of the 
fact that other very popular algebraic manipulations systems are on 
the market (like Mathematica or MAPLE) thus the area of people in- 
terested in algebraic programming routines for Dirac equation should 
be much larger. In this perspective we developed similar programs and 



routines for MAPLE |^ platform using the package GRTensorll 



(adapted for doing calculations in General Relativity). Because there 
is no portability between the two systems we were forced to compose 
completely new routines, in fact using the same strategy we used in 
REDUCE : first, the Pauli and Dirac matrices algebra (using only the 
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MAPLE environment) and then the construction of the Dirac equation 
on curved spacetimes were the capabihties of GRTensorll package is 
used. Because the authors of GRTensorll offer also package versions 
for Mathematica we are sure that our MAPLE routines can be eas- 
ily adapted for Mathematica with the result of highly increasing the 
number of users of our product. 

This chapter is organized as follows : the next section presents a 
short review of the theory of Dirac fields on curved spacetime, point- 
ing out the main notations and conventions we shall use. The next 
section is devoted to a short overview of our routines and programs in 
REDUCE+EXCALC previously described in great detail in This 
section is necessary in view of the fact that we used the same strat- 
egy for constructing our programs in MAPLE as we pointed out above. 
Then the section nr. 4 contains a complete description of our programs 
in MAPLE+GRTensorIL We also included here some facts about the 
main differences (advantages and disadvantages) between the two al- 
gebraic programming platforms (REDUCE and MAPLE). Section 5 is 
devoted to the problem of including spacetimes with torsion in order 
to compute Dirac equation using our MAPLE procedures. The last 
section of the article includes a list of some of spacetimes examples we 
used in order to calculate the Dirac equation. Two of these examples 
are spacetimes with torsion thus we pointed out the contribution of 
torsion components to the Dirac field. Several applications of our pro- 
grams (in REDUCE or in MAPLE) are in view of our future research 
: searching for inertial effects in noninertial systems of reference (par- 
tially presented in |]33| for a Schwarzschild metric without torsion) or 
quantum effects (as in [Q) in order to provide new theoretical results 
for experimental gravity 



3.2 Pauli and Dirac matrices algebra and 
Dirac equation on curved space-time 

The main problem is to solve algebraic expressions involving the Dirac 
matrices |^ , |]2^ , ||2^ . To this end it is convenient to construct explic- 
itly these matrices as a direct product of several pairs among the Pauli 
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matrices ai,i = 1, 2, 3, and the 2x2 unit matrix. Thus all the calcula- 
tions will be expressed in terms of the Pauli matrices and 2-dimensional 
Pauli spinors. Consequently the result will be obtained in a form which 
is suitable for physical interpretations. We shall consider the Pauli ma- 
trices as abstract objects with specific multiplication rules. Thus we 
work with operators instead of their matrices in a spinor representa- 
tion. However, if one desires to see the result in the standard Dirac 
form with 7 matrices it will be sufficient to use a simple reconstruction 
procedure which will be presented in a next section. 

We shall consider the Dirac formalism in the chiral form where the 
Dirac matrices are 12^ : 



7 



The Dirac spinor : 



" 


1 ' 




" a,' 




"-10" 


1 





_-ai 


, 7^ = 


1 



e n 



D 



(3.1) 



(3.2) 



involves the Pauli spinors ipi and which transform according to the 
irreducible representations (1/2,0) and (0,1/2) of the group SL(2,C). 
In this representation the left and right-handed Dirac spinors are 



1 -7^ 







1 + 7^ 



R 







(3.3) 



and, therefore, the Pauli spinors (pi and will be the left and the 
right-handed parts of the Dirac spinor. The SL(2,C) generators are 



(3.4) 



It is clear that 1-Ld = 'H ® IH. { where 7i is the two-dimensional 
space of Pauli spinors) and, therefore the Dirac spinor can be written 
as: 



1 




and ^2 




1 



(3.5) 



^' = 6 ® V^/ + 6 ® V^r with ^1 : 
while the 7-matrices and the SL(2,C) generators can be put in the form 

1 (3.6) 
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(3.7) 



These properties allow to reduce the Dirac algebra to that of the Pauli 
matrices. 

In order to introduce Dirac equation on curved space-time we have 
always used an anholonomic orthonormal frame because at any point 
of spacetime we need an orthonormal reference frame to describe the 
spinor field as is already pointed before |^). The Dirac equation in 
a general reference of frame, defined by an anholonomic tetrad field is 



and where S^'^ are the SL(2,C) generators ( p. 41) and F^p^ are the an- 
holonomic components of the connection. 



3.3 Review of the REDUCE+EXCALC 
routines for calculating the Dirac equa- 



We shall describe here those part of the program realizing the Pauli 
and Dirac matrix algebras. In the first lines of this sequence we in- 
troduce the operators and the non-commuting operators being useful 
throughout the entire program. The Pauli matrices are represented us- 
ing the operator p with one argument. The Dirac matrices are denoted 
by gam of one argument (an operator if we use only REDUCE, or for 
EXCALC package it will be a 0-form with one index) while the oper- 
ator dirac stands for the Dirac equation. The SL(2,C) generators are 
denoted by the 0-form s(a,b). The basic algebraic operation, the com- 
mutator (com) and anticommutator (acom) are then defined here only 
for commuting (or anticommuting) only simple objects ("kernels"). For 
commuting more complex expressions, (in order to introduce some nec- 
essary commutation relations) a more complex operator is necessary to 





(3.9) 



tion 
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introduce. Other objects, having a more or less local utilization in the 
program will be introduced with declarations and statements at their 
specific appearance. 

The main part of the program is the Pauli subroutine : 

LET p(0)=l; 

LET p(2)*p(l)=-p(l)*p(2) ; 
LET p(l)*p(2)=i*p(3); 
LET p(3)*p(l)=-p(l)*p(3); 
LET p(l)*p(3)=-i*p(2); 
LET p(3)*p(2)=-p(2)*p(3); 
LET p(2)*p(3)=i*p(l) ; 
LET p(l)**2=l; 
LET p(2)**2=l; 
LET p(3)**2=l; 

The Pauli matrices, (Tj appear as p(i) while the 2x2 unity matrix is 
p(0)=l. The properties of the Pauli matrices arc given by the above 
sequence of 10 lines. The direct product denoted by pd operator has 
the properties introduced as : 

for all a,b,c,u let pd(a,b)*pd(c,u)=pd(a*c,b*u) ; 

for all a,b let pd(a,b)**2=pd(a**2,b**2) ; 

for all a,b let pd(-a,b)=-pd(a,b) ; 

for all a,b let pd(a,-b)=-pd(a,b) ; 

for all a,b let pd(i*a,b)=i*pd(a,b) ; 

for all a,b let pd(a, i*b)=i*pd(a,b) ; 

for all a let pd(0,a)=0; 

for all a let pd(a,0)=0; 

let pd(l,l)=l; 

Some difficulties arise from the bilinearity of the direct product which 
requires to identify all the scalars involved in the current calculations. 
This can be done only by using complicated procedures or special as- 
signments. For this reason we shall use a special definition of the direct 
product (pd) which gives up the general bilinearity property. The op- 
erator pd will depend on two Pauli matrices or on the Pauli matrices 
with factors — 1 or ±i. It is able to recognize only these numbers but 
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this is enough since the multiphcation of the Pauh matrices has just 
the structure constants ±1 and ±i (we have a^aj = 6ij + ieijk<Jk)- 

Thus by introducing the multiphcation rules of the direct product it 
will be sufficient to give some instructions (see above) which represent 
the bilinearity defined only for the scalars —1 and i. The next two 
instructions represent the definition of "0" in the direct product space 
while the last one from the above sequence introduces the 4x4 unit 
matrix. 

The 7" matrices can be defined now with the help of our direct 
product; also we added here the definition of the SL(2,C) generators 
from eq. ( P^ ) : 

gamd) :=i*pd(p(2) ,p(l)) ; 
gaiii(2) :=i*pd(p(2) ,p(2)) ; 

gam(3) :=i*pd(p(2) ,p(3) ) ; % Remember that gam's are 

gaiii(0):= pd(p(l),l); % 0-forms !!! 

gam5 := -pd(p(3) ,1) ; 7. instead of "gam(5) " 



s(a,b):= i*com(gaiii(a) ,gam(b) )/4; 

All the above program lines we have presented here can be used for 
the Dirac theory on the Minkowski space-time in an inertial system of 
reference. Here we shall first point the main differences which appear 
when one wants to run our procedures on some curved space-times or 
in a noninertial reference of frame. Some of these minor differences are 
already integrated in the lines presented in the precedent section. 

First of all we must add, at the beginning of the program some 
EXCALC lines containing the metric statement. We have always used 
an anholonomic orthonormal frame because at any point of spacetime 
we need an orthonormal reference frame in order to describe the spinor 
field as is already pointed before ||3^). 



After the above metric statement we must add in the program all 
the procedures described in the last section. Then we can introduce 
some lines calculating the Dirac equation in this context. As a result 
we have used the next sequence of EXCALC lines : 



pform {der(j) ,psi}=0; fdomain psi=psi(x,y,z,t) ; 
der(-j):= ee(-j)_|d psi + (i/2)*s(b,h)*cris(-j ,-b,-h) ; 
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operator derpO , derpl , derp2 , derp3 ; 
noncom derpO , derp 1 , derp2 , derp3 ; 

let @(psi ,t)=derpO ; 
let @(psi ,x)=derpl ; 
let (§(psi ,y)=derp2 ; 
let @(psi,z)=derp3; 

dirac := i*has*gain(j)*der(-j)-m*c; 

ham:= -(gain(0)*(l/(ee(-0)_ I d t) ) *dirac-i*has*derpO) ; 

In defining the Dirac derivative der we have introduced also an formal 
Dirac spinor (psi) being a 0-form and depending on the variables im- 
posed by the symmetry of the problem. It is just an intermediate step 
(in fact a "trick") in order to obtain the partial derivative components 
as operators, because after calculating the components of the covariant 
derivative (der(-j)- see above) we have to replace the partial derivatives 
of psi with four non-commuting operators derpO, derpl ... derp3. 
The Dirac operator is thus defined as dirac:= {ih^^D^ — mc)^ and 
finally the Dirac Hamiltonian (ham) is obtained from the canonical 
form of the Dirac equation : 

at 

which we shall use later, in the study of the nonrelativistic approxima- 
tion of the Dirac equation in noninertial reference frames. 

The results we have obtained after processing the program lines we 
have presented until now contains only the Pauli matrices and direct 
products of Pauli matrices. When one wish to have the final result in 
terms of the 7-matrices and SL(2,C) generators (and not in terms of 
direct products of Pauli matrices) the procedure rec will be used : 

operator gama,gen; 
noncom gama,gen; 
PROCEDURE rec (a); 
begin; 
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wsl :=sub(pd(p(l) , l)=gama(0) ,a) ; 
wsl :=sub(pd(p(2) , l)=-i*gama(0) *gama(5) ,wsl) ; 
""1 :=sub(pd(p(3) , l)=-gama(5) ,wsl) ; 
l:=sub(pd(l,p(l))=2*gen(2,3) ,wsl) 
l:=sub(pd(l,p(2))=2*gen(3,l) ,wsl) 
l:=sub(pd(l,p(3))=2*gen(l,2) ,wsl) 



ws 
ws 
ws 
ws 

for k:=l:3 do 



«wsl :=sub(pd(p(2) ,p(k) )=-i*gama(k) ,wsl) ; 

wsl : =sub(pd(p(l) ,p(k) )=gama(k) *gama(5) , wsl) ; 
wsl :=sub(pd(p(3) ,p(k) )=2*i*gen(0,k) ,wsl)»; 
return wsl; 
end; 

This is an operator depending on an expression involving matrices (a) 
which reconstructs the 7-matrices and the SL(2,C) generators from the 



direct products of Pauh matrices according to eq. ^]6| and pTf . 

As a very important remark we must point out that the new intro- 
duced operators gama and gen does not represent a complete algebra. 
They are introduced in order to have the result in a comprehensible 
form. Thus this form of the result cannot be used in further compu- 
tations. Only the results obtained before processing the rec procedure 
can be used, in order to benefit of the complete Pauli and Dirac ma- 
trices algebra. We have used this rec procedure only for pointing out 
our results in a more comprehensible form. 



3.4 MAPLE+GRTensorll procedures for 
calculating the Dirac equation 

Here we shall present, in details, our procedures in MAPLE-|-GRTensorII 
for calculating the Dirac equations, pointing out the main differences 
between MAPLE and REDUCE programming in obtaining the same 
results. The first major problem appears in MAPLE when one try to 
introduce the Pauli and Dirac matrices algebra. In MAPLE this will 
be a difficult task because the ordinary product (assigned in MAPLE 
with "*") of operators is automatically commutative, associative, lin- 
ear, etc. like an ordinary scalar product - in REDUCE these properties 
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are active only if the operators are declared previously as having such 
properties. Thus we have to define two special product operators : for 
Pauli matrices (Tq,q; = 0, 1,2,3 (assigned in our procedures with pr) 
and for the direct product of Pauli matrices (assigned here also with 
the operator pd(",")) which is assigned with "&;p". As a consequence 
we have to introduce long lists with their properties as : 

> def ine(sigma,sigma(0)=l) ; 

> def ine(pr ,pr (1 , 1)=1 ,pr (1 , sigma(l) )=sigma(l) ,pr (1 , sigma(2) )= 
sigma(2) ,pr (1 , sigmaO) )=sigma(3) ,pr (sigma(l) , l)=sigma(l) , 
pr(sigma(2) , l)=sigma(2) ,pr(sigma(3) , l)=sigma(3) , . . . 
pr(sigma(l) ,sigma(l))=l,pr(sigma(2) , sigma(2) )=1 , . . . 

> def ine(pd,pd(0,a: :algebraic)=0,pd(a: :algebraic,0)=0,pd(l, 1)=1, 
pd(I*a: : algebraic, b: : algebraic) =l*pd (a, b) , 

pd(-I*a: :algebraic,b: : algebraic) =-I*pd (a, b) , 
pd(a: : algebraic , I*b : : algebraic)=I*pd(a,b) , 
pd(a: : algebraic , -l*b : : algebraic)=-l*pd(a,b) , 
pd(-a: : algebraic, b: : algebraic) =-pd (a, b) , 
pd(a: : algebraic, -b: : algebraic) =-pd (a, b) ) ; 

> def ineC&p' , '&p' (-a: :algebraic,-b: : algebraic) = '&p' (a,b) , . . . 

We dropped the complete hst of the properties of the special products 
pr and 8zp, being very long. Of course the reader may ask why is not 
much simpler to declare, as an example the &j9 as being linear (or mul- 
tilinear) ? Because in this case the operator does not act properly, the 
linearity property picking out from the operator all the terms, being 
or not Pauh matrices or direct products pd of Pauh matrices. Thus is 
necessary to forget the linearity and to introduce, as separate properties 
all the possible situations to appear in the calculus. The result is that 
the program become very large with a corresponding waste of RAM 
memory and speed of running. This will be the main disadvantage of 
MAPLE version of our program in comparison with the short (and, 
why not, elegant) REDUCE procedures. Of course, in a more compact 
version of our programs, we defined MAPLE routines with these op- 
erators, and the user need only to load at the beginning of MAPLE 
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session these routines, but there is no significative economy of memory 
and running time. 

The next step is to define Dirac 7-matrices and a special commuta- 
tor (with &p) : 

> define (gam, 

gam(l)=I*pd(sigma(2) ,sigma(l)) , 
gam(2)=I*pd(sigma(2) ,sigma(2)) , 
gam(3)=I*pd(sigma(2) ,sigma(3)) , 
gam(0)=pd(sigma(l) ,1) ,gam(5)=-pd(sigma(3) ,1)) ; 

> def ine(comu,comu(a: :algebraic,b: :algebraic)=a &p b - b &p a); 

The next program-hnes are in GRTensorll environment. For this is 
necessary to load previously the GRTensorll package and then to load 
the corresponding metric (with qload(...) command. It follows then : 

> grdef('SS{ "a ~b }'); 

> grcalc (SS (up , up) ) : 

> (I/4)*comu(gam(0) ,gam(0)) ; 

> (I/4)*comu(gam(l) ,gaiii(0)) ; 



> (I/4)*comu(gam(3) ,gaiii(3)) ; 

> grdisplay(SS(up,up)) ; 

> grcalc (Chr (dn, dn, dn) ) ;grdisplay(Chr(dn,dn,dn)) ; 

> grcalc (Chr (bdn,bdn,bdn) ) ; grdi splay (Chr (bdn,bdn,bdn) ) ; 

These are a sequence of commands in GRTensorll for defining the 
SL(2,C) generators S^^ (as the tensor SS{ "a "b}) using formula ( |3.4| ) 
and for the calculus of Christoffel symbols in an orthonormal frame 
base (Chr(bdn,bdn)). Here is active one of the main advantages of 
MAPLE+GRTensorll platform, namely the possibility of the calculus 
of the tensor components both in a general reference frame or in an 
anholonomic orthonormal frame which is vital for our purpose of con- 
struction the Dirac equation. 
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Next we have to define, as two vectors the Dirac-7 matrices (as- 
signed as the contravariant vector ga{ "a } and the derivatives of the 
wave function ip (assigned as the covariant vector Psid{ a } in order 
to use the facihties of GRTensorll to manipulate with indices : 

> grdef('ga{ "a }: = [gain(0) ,gajii(l) ,gam(2) ,gam(3)] ') ; 

> grdisplay(ga(up)) ; 

> grdef('Psid{ a }:=[diff (psi(x,t) ,t) ,diff (psi(x,t) ,x) , 

diff (psi(x,t) ,y) ,diff (psi(x,t) ,z)] ') ; 

> grcalc(Psid(dn) ) ;grdisplay(Psid(dn) ) ; 

> grcalc(Psid(bdn) ) jgrdisplay (Psid(bdn) ) ; 

The next step defines a term which will be the term i^^^'^Vyp^ from 
equation ^]9| : 

> grdef('de{ i }:=(I/2)*SS{ "a "b }*Chr{ (i) (a) (b) }'); 

> grcalc(de(dn) ) ;grdisplay(de(dn) ) ; 

Observing that the components of de{ i } are polynoms containing 
direct products pd(...) of Pauli matrices and the fact that the prod- 
uct between Y and i^S^^Vyp^ from equation ^]8| is, in fact the special 
product &p we have to obtain the term SP^V^p^ (denoted below with 
the operator dd) by a special MAPLE sequence which in fact split 
the components of de{ i } in monomial terms and then execute the 
corresponding &p product, finally reconstructing the dd operator : 

aO : =expand (gr component (de (dn) , [t] ) ) ; aOO : =0 ; 

> uO:=whattype(aO) ;uO; nops(aO); 

> if uO='+' then for i from 1 to nops(aO) do 

aOO : =aOO+I*h*grcomponent (ga(up) , [t] ) &p op(i,aO) od else 
aOO : =I*h*grcomponent (ga(up) , [t] ) &p aO fi; aOO; 

> al : =expand (gr component (de (dn) , [x] ) ) ; al 1 : =0 ; 

> ul : =whattype(al) ;ul ; 

> nops (al) ; 

> if ul='+' then for i from 1 to nops(al) do 

all : =all+l*h*grcomponent (ga(up) , [x] ) &p op(i,al) od else 
all : =I*h*grcomponent (ga(up) , [x] ) &p al fi;all; 
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> grdef Cdd') ; grcalc(dd); 

> a00+all+a22+a33; 

> grdi splay (dd) ; 

Finally the Dirac equation is obtained as : 

> grdef ('dirac := 

I*h*ga{ "1 }*Psid{ (1) } + dd*psi(x,t) - in*c*psi(x,t) ') ; 

> grcalc (dirac) ; 

> grdisplay (dirac) ; 

In order to obtain the Dirac equation in a more comprehensible form 
we have the next sequence of MAPLE commands (similar to the recon- 
struction rec procedure from the REDUCE program : 

> define ( 'gen' ) ; 

> define (' gama' ) ; 

> grmap (dirac, subs, pd(sigma(l) , l)=gaiiia(0) ,pd(sigma(2) , 1)=-I*gaina(0) * 
gama (5) ,pd(sigma(3) , l)=-gama(5) ,pd(l , sigma(l) )=2*gen(2 ,3) ,pd(l , 
sigina(2))=2*gen(3,l) ,pd(l,sigma(3))=2*gen(l,2) ,pd(sigma(2) ,sigma(l))= 
-I*gaiiia(l) ,pd(sigma(2) , sigma(2) )=-I*gama(2) ,pd(sigma(2) ,sigina(3)) = 
-I*gaiiia(3) ,pd(sigma(i) ,sigma(i))=gama(i)*gama(5) ,pd(sigma(i) ,sigma(2)) = 
gama(2) *gaina(5) ,pd(sigma(l) , sigma(3) )=gama(3) *gama(5) ,pd(sigma(3) , 
sigma(l))=2*I*gen(0,l) ,pd(sigina(3) ,sigina(2))=2*I*gen(0,2) ,pd(sigma(3) , 
sigina(3))=2*I*gen(0,3) , 'x') ; 

> grdisplay (dirac) ; 

where, of course, as in REDUCE version, the operators gen and 
gama does not represent a complete algebra. 

3.5 Dirac equation on spacetimes with 
torsion and computer algebra 

We shall present here how we adapted the already presented in the pre- 
vious section our MAPLE+GRTensorll programs in order to calculate 
the Dirac equation on space-times with torsion. 
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The geometrical frame for General Relativity is a Riemannian space- 
time but one very promising generalization is the Riemann-Cartan ge- 
ometry which (i) is the most natural generalization of a Riemannian 
geometry by allowing a non-symmetric metric-compatible connection, 
(ii) treats spin on the same level as mass as it is indicated by the group 
theoretical analysis of the Poincare group, and (iii) arises in most gauge 
theoretical approaches to General Relativity, as e.g. in the Poincare- 
gauge theory or supergravity However, till now there is no 

experimental evidence for torsion. On the other hand, from the lack of 
effects which may be due to torsion one can calculate estimates on the 



maximal strength of the torsion fields In this aspect we think that 
is possible, using computer algebra facilities to approach new theoret- 
ical aspects on matter fields (for example the Dirac field) behavior on 
spacetimes with torsion in order to point out new gravitational effects 
and experiments at microscopic level. 

A metric compatible connection components in a Riemann-Cartan 
theory is related to the torsion components by (see - eq. (1.18)) 

(3.10) 

where T^p^^ are the components of the riemannian connection, 
is the object of anholonomicity and Tap'y are the components of the 
torsion. The idea is to replace the connection components from the 
covariant derivative appearing in eqs. ( |3.8[ - p.9| ) with the above ones, of 
course after calculating then in an orthonormal anholonomic reference 
frame suitable for calculation of the Dirac equation. Thus the sequence 
for calculating the de{dn} operator (see above) should be replaced 
with 

> grdef('de{ i }:=(I/2)*SS{ 'a 'b }*(CHR{ (i) (a) (b) }+ 

(l/2)*(tor{ (i) (a) (b) } -tor{ (a) (b) (i) } 
+ tor{ (b) (i) (a) })))'); 

> grcalc (de (dn) ) ; grdisplay(de(dn) ) ; 

where the new connection components CHRa, b, c are now defined 
by the sequence 



> grdef('ee{ a ~b }:= wl{ "b >*kdelta{ a $x } 
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+ w2{ "b >*kdelta{ a $y > + 
w3{ ~b }*kdelta{ a $z > + w4{ "b }*kdelta{ a $t }'); 

> grcalc(ee(dn,up) ) ; 

> grdisplay(ee(dn,up)) ; 

> grdef('CC{ a b c }:=2*ee{ a "i >*ee{ b " j }*ee{ c [j ,i] >'); 

> grcalc(CC(dn,dn,dn) ) ; grdisplay(CC(dn,dn,dn) ) ; 

> grdef('CHR{ (a) (b) (c) > := Chr{ (a) (b) (c) > - 

(1/2)*(CC{ a b c > - 
CC{ b c a } + CC{c a b }) ') ; 

> grcalc(CHR(bdn,bdn,bdn)) ; grdisplay(CHR(bdn,bdn,bdn)) ; 

and the rest of the routines are unchanged. The only problem re- 
mains now to introduce in an adequate way the components of the 
torsion tensor. We used a suggestion from pointing that we can 
assume that the 2-form T" associated to the torsion tensor should have 
the same pattern as the c/^^'s where 6'" is the orthonormal coframe, 
who's components are denoted in GRTensorll with wl{dn} ...w4{dn}. 
Thus this operation it is possible only after we introduced the metric 
(with qload command). Calculating then the derivatives of the or- 
thonormal frame vector basis components we can introduce the torsion 
components by inspecting carefully these derivatives. Here is an ex- 
ample of how this is possible in MAPLE+GRTensorll using one of the 
metric examples presented in the next section: 

> grcalc(wl (dn,pdn) ) ;grcalc(w2(dn,pdn) ) ; 

> grcalc(w3(dn,pdn) ) ;grcalc(w4(dn,pdn) ) ; 

> grdisplay (wl (dn,pdn) ) ; grdisplay(w2(dn,pdn) ) ; 

> grdisplay (w3 (dn, pdn) ) ; grdisplay (w4(dn,pdn) ) ; 

> grcalc (wl (bdn , pbdn) ) ; grcalc (w2 (bdn , pbdn) ) ; 

> grcalc (w3 (bdn, pbdn) ) ; grcalc (w4 (bdn, pbdn) ) ; 

> grdisplay (wl (bdn, pbdn) ) ; grdisplay (w2 (bdn, pbdn) ) ; 

> grdisplay (w3 (bdn, pbdn) ) ; grdisplay (w4 (bdn, pbdn) ) ; 

> grdef('tor-[ "a b c }:=wl{ b ,c }*kdelta-[ ~a $x } + 

w2{ b ,c }*kdelta{ ~a $y } + 
w3{ b ,c }*kdelta{ "a $z } + 
w4{ b ,c }*kdelta{ "a $t }'); 

> grcalc (tor (up , dn , dn) ) ; grdisplay (tor (up , dn , dn) ) ; 

> grmap(tor (up,dn,dn) ,subs,f (x)=v4(x) ,h(x)=v3(x) ,g(x)=v2(x) , 'x') ; 
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> grdisplayCtor (up,dn,dn) ) ; 

> grcalc(tor(bup,bdn,bdn) ) ; 

> grdisplay(tor(bup,bdn,bdn)) ; 

The reader can observe that we first assigned the components of the 
torsion tensor (here denoted with tor{up,dn,dn}) then after display- 
ing his components we ca decide to substitute new functions describing 
the torsion instead of the functions describing the metric. Of course 
finally we calculate the components of the torsion in an orthonormal 
anholonomic reference frame (tor{bup,bdn,bd}). 



3.6 New results 

This section is devoted to a list of more recent results we obtained by 
running our procedures in MAPLE+GRTensorll already described in 
the previous sections. First we tested our programs by re-obtaining the 
form of the Dirac equations in several spacetime metrics we already 
obtained with REDUCE+EXCALC procedures and reported in our 



articles The concordance of these results with the previous 

ones was a good sign for us to proceed with more complicated and new 
examples, including ones with torsion. Here we shall present some of 
these examples. 

1. Conformally static metric with $ and S constant where the line 
element is : 

Thus the Dirac equation becomes : 

the ^ 7M — ^— - ^ , , b'{r) - — — + 7 + 



a{r) dt 2a{r)b{r) a{r)r J \r dr r sin{9) dO 

1 . 1 n/3 d 

ip — mcip = U 



where h'{r) is the derivative db{r)/dr. 
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2. Taub-NUT spacetime having the hne element as : 
ds'^ = -Al^U{t)dy^ - 8l^U{t) cos{e)dyd(j) - {t^ + f)de'^ 

+ {-AfU{t) cos{ef - {t" + /2) sin{ef)d(t)'' + j^dt' 

U [t) 

the coordinates being {y,9,phi,t). We obtained the Dirac equation as 

ih 



t^ + P 



-7 



\A^/TJ{t) y ' ' \ dt 



Of ^ U'(t) + Ju{t)[l + ^\\-cotg{e)^Jl^^'' 



m 

mcil){t) = 



3. Gddel spacetime, having the hne element as : 

ds^ = -a^dx^ + ia^e^^c??/^ + 2a^e'^cdydt - a^dz^ + a?c^dt^ 

in coordinates {x,y,z,t). Here the Dirac equation is simply : 
1 



ih 
a 



\2 dx J ^ Qy cdtj ^ dz cdt 



i){x,y,z,t) 
—mcip{x, y, z,t) =0 

4. McCrea static spacetime with torsion having the line element 
as 

in (x, y, 2;, t) coordinates. If the coordinate lines of y are closed with < 
y < 2tt and — oo<2<oo,0<x<oo, the spacetime is cylindrically 
with y as the angular, x the cylindrical radial and z the longitudinal 
coordinate. If — 00, x,y,z < 00 the symmetry is pseudo-planar. In 
26| McCrea considers the simplest solution of Einstein equations with 



cosmological constant as 

f = h = h = qx/3 (3.11) 

and the cosmological constant turns to be We shall first consider 

the general case specializing the results at the final step of the program 
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to the above particular solution. Running our MAPLE+GRTensorll 
procedures first we obtain the torsion tensor component as : 



dx ' dx ' dx 



the rest of the components being zero. This time we have obtained the 
Dirac equation, depending also on the components of the torsion tensor 
as : 

•^1^1 fv2(x) dv2ix) ^ ^^3(^) dv3{x) ^ ^M(^) gf4(x) _ dfjx) _ dgjx) 
2 \ dx dx dx dx dx 

dh{x) ^ 2^--] iIj{x) - mcipix) = 



dx dx ^ 

Of course when we specialize to the particular solution proposed by 
McCrea in we have to assign the form of metric functions as in 
and we can then take the torsion functions as 

v2 = v3 = vi = v{x) 

and the Dirac equation becomes 

5 Schwarzschild metric with torsion. This example it is interesting 
in view of recent investigations on the contribution of torsion in grav- 
ity experiments using atomic interferometry [^. Here we used the 
Schwarzschild metric having the line element written as 

ds^ = e^^'Ur^' + r^de^ + sin{efd<p' - e'^^'^dt' 

Here we preferred to specialize the form of A(r) and z/(r) functions as 

2m 1 
v[r) = 1 



r A(r) 

(here we have G = c = 1) later after obtaining the form of the Dirac 
equation in term of A(r) and z/(r) functions. 
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Using the same "trick" as in the previous example, we choose the 
components of the torsion tensor as 

T^r = ^^eV^/UO ; Tl = 1 ; T^ = s^n{e) 



= ^^eV^/^« ; Tte^rcosie) 

the rest of the components being zero. Running away our procedures 
we obtain the Dirac equation containing terms with torsion tensor com- 
ponents as : 



ih 



-A(r)/2^1 ( '9/2 (r 



7 



,/2(r-)/2 _ 



y dr 



dv{r) 
dr 



4 d\ 
- + 2(l + sm(e))+4— + 
r or 



--f{cos{e)-cotg{e)) 



■0(r) — mip{r) — 
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